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Abstract 

Neutrinos from far-away sources annihilating at the Z resonance 
on rehc neutrinos may give origin to the extreme-energy cosmic rays 
(EECR). If "Z-bursts" are responsible for the EECR events, then we 
show that the non-observation of cosmic ray events at energies above 
2 xlO^'' eV by the AGASA Collaboration implies a lower bound 
~ 0.3 eV on the relic neutrino mass. Since this mass exceeds the 
mass-squared differences inferred from oscillation physics, the bound 
in fact applies to all three neutrino masses. Together with the upper 
bound provided by comparisons of the CMB anisotropy with large- 
scale structure, this bound leaves only a small interval for neutrino 
masses around 0.3 eV, if Z-bursts are to explain the existing EECR 
events. 
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1 Introduction 



The existence of extreme-energy cosmic rays (EECR) with energies above 
the Greisen-Zatsepin-Kuzmin (GZK) cutoff [1] of about 5 x 10^^ eV, presents 
an outstanding problem [2]. Nucleons and photons with those energies have 
short attenuation lengths and could only come from distances of 100 Mpc or 
less [3, 4], while plausible astrophysical sources for those energetic particles 
are much farther away. 

Data from the HiRes experiment [5] brought the violation of the GZK 
cutoff into question. Yet the EE events in even the HiRes data set remain 
unexplained since the local Universe (~ 100 Mpc) is devoid of strong can- 
didate sources. This controversy will be solved conclusively by the Pierre 
Auger hybrid observatory [6], perhaps as soon as summer of 2005. Here 
we assume that the published AGASA spectrum is correct. Interestingly, it 
is the absence of events above 2 x 10^° eV in these data, rather than the 
presence of events above Eqzk ~ 5 x 10^^ eV, that motivates this work. 

Among the solutions proposed for the origin of the highest energy events 
observed by AGASA, an elegant and economical solution to this problem is 
the "Z-burst" mechanism: annihilation at the Z-resonance of extreme-energy 
neutrinos (z/ee) coming from remote sources, and relic background neutrinos 
within ~ 50 Mpc of Earth, produces the nucleon and photon EECRs [7]. 
The observed EECRs from the Z-bursts are the emission analogues of the 
absorption features (Z-dips) predicted long ago [8]. One of the most appealing 
features of the Z-burst mechanism is that the energy scale of ^ lO^'' eV at 
which the unexpected events have been detected, is generated naturally. The 
Z-resonance occurs when the energy of the incoming iz-^e is 



where m„ is any of the three masses of the relic neutrinos. Given the lower 
limit ~ 0.04 eV deduced from atmospheric oscillations for the heaviest neu- 
trino mass, at least one E'res is below 10^^ eV. 

Since the individual energies of the nucleons and photons emerging from 
the Z-burst cannot exceed the total energy of the burst, £^res is the new 
end-point of the EECR energy-spectrum in this mechanism. Partitioning 
this burst energy among the (N) ~ 40 final state particles, one arrives at 
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precisely the primary energies needed to produce events observed above the 
GZK energy. 

Combining their recent measurements of the anisotropies of the Cosmic 
Microwave Background (CMB) radiation with data on the large scale struc- 
ture of the Universe, WMAP [9] produced a strong limit on neutrino masses, 
T^irriy^ < 0.69 cV. Since a single neutrino mass above ~ 0.04 cV implies near 
mass-degeneracy for all three active neutrinos (given the mass-squared split- 
tings from neutrino oscillation data, described below) one has rrii, < 0.23 eV 
at the 95% CL. However, objections [10] to the priors assumed, or to the data 
sets included, have led to a relaxed bound Simjy^ < 1 eV, or < 0.33 eV. 
A subsequent analysis [11] which includes WMAP, 2dF, and SDSS data, and 
another which includes these and galaxy cluster data [12], have arrived at the 
bound, Ejm^,. ^ 0.7 eV. Still another analysis [13] with CMB and LSS data 
finds T^iirLy^ ^ 1 eV, but finds 'EiUiui ^ 0.6 eV when priors from supernova 
data and the Hubble Key Project are included. These newer bounds are very 
similar to the original WMAP bound. 

In this work we focus on a particular feature of the AGASA spectrum, 
namely the end-point energy. Requiring that the Z-burst mechanism not 
overproduce events beyond this AGASA end-point, we derive a lower bound 
on the neutrino mass. The significance of an end-point energy for constrain- 
ing model fluxes has been noted much earlier in the context of topological 
defect decay [14]. A more ambitious project would be to actually fit the 
neutrino mass to the full AGASA spectrum. This has been done, by Fodor, 
Katz and Ringwald (FKR) [15]. However, the allowed range of the neutrino 
mass that results is not tight. Rather, it depends sensitively on how one pa- 
rameterizes the transition from non-burst spectrum to burst spectrum near 
and above the ankle. FKR find ±2a fitted mass ranges of 0.1 to 1.3 eV if 
all EE events are assumed to originate with Z-bursts, 0.02 to 0.8 eV if an 
additional extra-galactic source of EE protons is allowed, and 1 to 7 eV if an 
additional Galactic-halo source of EE protons is allowed. 

Super-Kamiokande has provided a strong evidence for the oscillation in 
atmospheric showers of two neutrino species with mass-squared splitting Sm^ 
— m\ — m\ = (1 — 3) X 10~^ eV, consisting of nearly equal amounts of and 
Vr-i and little or no [16]. If neutrino masses are hierarchical, like the other 
leptons and quarks, then \/5'w? ~ 0.04 eV = msK is the mass of the heavier 
of the two oscillating neutrinos (call it vsk)- Some previous works [17, 18] on 
Z-burst production of the EECRs have made this assumption. However we 
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will show in this work that ~ 0.04 cV is not compatible with the AGASA 
end-point spectrum if Z-bursts produce the EECR. The reason is that with 
nil, — 0.04 eV, the new EECR cutoff E^es — 10^^ eV predicts too many super 
GZK EECR events beyond the AGASA end-point. 

As the relic neutrino mass increases, Ej-^s oc 1 / decreases, and the fea- 
tures in the EE spectrum move progressively to lower energies. In particular, 
the number of events predicted at high energy decreases, and the number be- 
yond the AGASA end-point becomes compatible with zero for rrii, ^ 0.3 eV. 
Note that our result is much sharper than, but compatible with, the allowed 
mass ranges of the FKR models with no additional EE proton source, or with 
an additional extra-galactic source of EE protons. 

In the next section we present our simulations and the resulting spectrum 
of EECR in detail. 

2 Spectrum of EECR from Z-bursts 

We have performed simulations of the photon, nucleon and neutrino fluxes 
coming from a uniform distribution of "Z-bursts", namely annihilations 
at the Z pole {I'l' — > Z — > PI---)- The burst energy is given in Eq. 1. Our 
simulations cover the range msK < rrii, < 1 eV for relic neutrino mass. The 
Z-bursts were simulated using PYTHIA 6.125 [19], and the absorption of 
photons and nucleons was modeled using encrgy-attcnuation lengths provided 
by Bhattacharjee and Sigl [20], supplemented by radio-background models 
by Protheroe and Biermann [21]. 

We simulated a uniform distribution of about 10^ Z particles up to a 
maximum Zmax — 2. If other Z-burst spatial distributions are used, the 
main features of the spectrum of the ultra-high energy nucleons, photons 
and neutrinos remain the same as given below. 

The decay of the Z bosons through all possible channels was done au- 
tomatically by the PYTHIA routines [19], using the default options of this 
program. For comparison with our later figures, we show in Fig. 1 the spec- 
tra given by PYTHIA, normalized per Z boson, for rrii, = 0.3 eV; redshifts 
and energy absorption arc not included in PYTHIA. The multiplicities that 
PYTHIA gives per Z-decay are 1.6 nucleons, 17 photons, 15 Ug, 30 and 
0.23 (in each case counting particles and antiparticles) . 

In our simulation, each Z boson generated by PYTHIA was placed on the 
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Figure 1: Spectra of stable particles produced per Z-decay by PYTHIA (no 
absorption or redshift included). The resonant energy in this example is 
Ej-es = 1.3 X 10^^ eV, coming from the choice rriy = 0.3 eV. 
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"event list" of the cascade generator at a randomly selected distance. The 
cascade of decay products was then boosted. The 7 factor corresponding to 
an energy E^, — E^^s in Eq. (1) is 



resulting from the boosted Z-decays. The gamma factor of the secondary 
particles was corrected to include their red-shifting subsequent to the Z- 
decay. 

Each neutrino, nucleon or photon was created by PYTHIA in the initial 
cascade at a fixed position, with fixed energy and direction of motion with 
respect to the Earth frame of reference. The distance each particle had to 
travel before reaching Earth was compared with the appropriate attenuation 
length in space for the particle energy. If the distance was smaller than the 
attenuation length, the particle was assumed to reach Earth unchanged. In 
the opposite case, the energy and momentum of the particle were degraded 
by a factor ^ after traveling a distance equal to the attenuation length (and 
the particle was placed again in the list constituting the cascade at the new 
position, with the degraded energy and momentum). 

Neutrinos do not interact in their propagation. Thus, the energy spectra 
of the three kinds of neutrinos were simply generated by summing the number 
of particles over energy bins and normalizing to neutrino multiplicities per 
Z times the total number of Z particles used. On the other hand, nucleons 
and photons do suffer energy- absorption interactions, and here we included 
energy absorption. The propagation process for nucleons and photons was 
continued until particles reached Earth and were counted in the final spectra, 
or until particle energies became too small to be significant. In the latter 
case, the particles were simply discarded from the cascade. At this point, 
the final nucleon and photon spectra were appropriately normalized to the 
total number of Z-particles used, as was done with the neutrinos. The results 
are given in Fig. 2 for rrii, = 0.3, with a fit to the six most energetic occupied 
bins in the AGASA data [22]. These bins consist of 24 events spanning the 
energy range from lO^^ *^ eV to 10^°'^ eV. 

The energy-attenuation of the nucleons and especially photons requires 
more discussion. For nucleon attenuation, uncertainties are small, and so the 





We then followed the propagation of the nucleons, photons and neutrinos 
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Figure 2: Predicted spectra for m,^ =0.3 eV from Z-bursts with a uniform 
distribution up to 2; = 2, added to a power law spectrum which fits the data 
below the ankle ~ 10^^ eV, and terminates at 4 x 10^^ eV. It is seen that 
EECR primaries above the ankle are nearly 100% nucleons up to lO^'' eV, and 
photons plus nucleons at higher energises. Also shown is the EE neutrino flux 
at the unique resonance energy which produces the required Z-burst rate. 
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Figure 3: As in Fig. 2 but with the flux multiphed by E^. The photon flux 
shown is the lowest of the three of Fig. 2 (labeled as (3)). 
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predicted GZK suppression is highly credible. For the nucleon energy at- 
tenuation, we used the length given by Bhattacharjee and Sigl in the Fig. 9 
of Ref. [20], based on results from Ref. [23] and [24]. However, for photon 
energy-attenuation, the length is quite uncertain, due to the uncertain spec- 
trum of the absorbing radio background. The photon flux shown in Fig. 2 as 
curve (1) results from using the attenuation length shown in Fig. 11 of Ref. 
[20], based on observations of Clark et al. [25]. Protheroe and Biermann [21] 
produced two models for the radio background which lead to shorter inter- 
action lengths (and therefore more absorption) than those based on Clark et 
al. From the interaction lengths they provided, we constructed approximate 
attenuation lengths for the models of Protheroe and Biermann. We did so 
by reducing the attenuation length based on Clark et al. by the ratio of the 
respective interaction lengths, and obtained the curves (2) and (3) in Fig. 2. 
Since mean interaction lengths and mean energy-attenuation lengths do not 
have exactly the same energy dependence, our construction is an approxi- 
mation. However, we believe the three curves which we display give a good 
representation of the possible range of predicted photon fluxes. 

In this work, we seek a bound on the excess of events predicted. There- 
fore, we adopt the lowest observable photon flux, curve number (3), associ- 
ated with the maximum photon absorption, when computing the total event 
number. This is the photon flux shown in Fig. 3. Notice that with this 
choice, the photon flux does not become comparable to the nucleon flux un- 
til ~ 2 x 10^° eV. We add to the Z-burst proton and photon fluxes a power 
law spectrum with slope —3.23, terminated somewhat arbitrarily (as it is not 
known where this component actually dies out) at 4 x 10^^ eV. Such a power 
law was found by AGASA to flt the data below the ankle, from 4 x 10^^ to 
10^9 eV (see Table V of [26]). 

3 Lower bound on the relic neutrino mass 

In Fig. 4 the observed and predicted integrated number of EECR event rates 
above 7 x 10^^ eV are shown for different values of the relic neutrino mass. 
The Z-burst fluxes have been normalized to reproduce the 24 highest-energy 
events observed by AGASA [22] in the energy range from 10^^'® eV to 
IQ20.4 gY These events are spread among the first six bins of Fig. 4. 

The area under each curve at the final bin shows the integrated number of 
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neutrino mass 


predicted event excess N^x 


Poisson probability e ^'''^ 


0.6 eV 


2.29 


10.1% 


0.5 eV 


2.76 


6.3% 


0.4 eV 


3.43 


3.2% 


0.3 eV 


4.43 


1.2% 


0.2 eV 


6.09 


0.23% 


0.1 eV 


9.67 


0.0063% 



Table 1: Poisson probabilities that no events are observed above the AGASA 
end-point, given the excess predicted by the Z-burst mechanism as a function 
of neutrino mass. 



events from 2.5 x lO^'^ eV to infinity predicted by that so- normalized Z-burst 
flux. This number of excess events can also be read off from the value of 
the curve at = 10^^'^ eV less the 24 normalizing events. Since AGASA 
has observed no events above 2.5x10^'^ eV although their aperture remains 
constant with increasing energy [22], it is justified to call the events predicted 
beyond 2.5x10^° eV an excess; we label the predicted number of excess events 
as A^ex- In Table 1 the predicted event excess above 2.5 x lO^'^ eV is shown 
as a function of the neutrino mass. Also shown is the Poisson probability 
that a fluctuation downward from the predicted value gives zero events, as 
observed. This latter quantity is just e~'^°^. 

The probabilities show that masses lighter than ~ 0.3 eV have a proba- 
bihty less than a percent. Expressed as Gaussian variances, one infers from 
the percentages in the table that ^ 0.2 eV is disfavored at ~ 3a while 
nil, > 0.4 eV is disfavored at only ~ 2a. Together with the cosmological 
limit on the neutrino mass discussed earlier, the value rrii, ~ 0.3 eV emerges 
as less than robust, but viable. 

We remark that our definition of "excess" as events above 2.5 x 10^° eV is 
somewhat conservative in that the two highest energy AGASA events occur 
at £; ~ 2 X 10^° eV (±25% experimental cncrgy-rcsohition), nearer to the 
beginning of the highest-energy occupied bin than to its end. On the other 
hand, we have not considered in our analysis the very-highest energy event, at 
3 X 10^° eV, from the Fly's Eye experiment. This is in the spirit of neglecting 
fiuorescent data, e.g. HiRes, as we have chosen to focus in a self-consistent 
way on just the ground-scintillator data set from AGASA. 
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Figure 4: Integrated Z-burst fluxes for different relic neutrino masses, nor- 
malized to reproduce the 24 highest-energy AGASA events in the energy 
range 10^^'^ to 10^°'^ eV (corresponding to the first six of seven bins shown 
here). The neutrino mass parameterizes the curves, with lighter mass at top 
and heavier mass at bottom. 1 1 



There are three immediate inferences to be drawn from the preferred 
neutrino mass value. The first is that since the energy of no secondary from 
the Z-burst may exceed the burst energy, the model predicts a new cutoff for 
EECRs. With 0.3 eV rehc neutrinos, the cutoff energy is -Ercs — 1.3 x 10^^ eV. 
Given the large mean-multiplicity in Z-decay, (N) ~ 40, a more realistic 
cutoff in the model is ~ 10^^ eV. Future experiments such as Auger, ANITA, 
RICE, EUSO, and SALSA may be able to test this prediction. The second 
inference is that with rrii, ~ 0.3 eV, a significant clustering of relic neutrinos 
around very large galactic super-clusters is expected [27]. This clustering 
could be evidenced as an enhancement of the EECR flux coming from the 
direction of M87, just 16 Mpc away, near the Virgo center. And the third 
inference is that given the smallness of mass-squared differences inferred from 
neutrino oscillation experiments, rrii, ~ 0.3 eV implies near mass-degeneracy 
for the three active neutrinos. Mass- degeneracy enhances the possibility of 
directly observing absorptive "Z-dips" in the EE neutrino spectrum [8]. 

4 Discussion 

We note that our analysis of the spectral shape of the AGASA data does not 
depend on the issue of rate. In particular, it does not depend on the EECR 
neutrino flux and it does not depend on the value of the relic neutrino target 
density. The latter could, in principle, be enhanced by either gravitational 
clustering or by a lepton asymmetry. In practice, simulations argue against 
gravitational clustering for light 1 eV) neutrinos on all but the very largest 
super-Galactic mass-scales [27] (which comprise a small sohd angle of sky), 
and BBN physics argues against a lepton asymmetry enhancement above 
~ 20% [28]. (Strictly speaking, the BBN limit is strongest for u^s, but the 
large lepton mixing angles inferred from oscillation studies extend the limit 
equally to i/^'s and i/^-'s.) Regardless of these kinds of rate arguments, our 
bound on the neutrino mass holds. 

Farrar and Piran [29] have argued that any mechanism accounting for the 
events beyond the GZK cutoff should also account for the events down to the 
ankle, including the observed isotropy and spectral smoothness. It is seen in 
our Fig. 2 that the Z-burst mechanism can accommodate the position of the 
ankle and all the events above it, if the position of the ankle is close to that 
measured by AGASA, at -Bankie = 10^^ ° eV (see [26], in particular Table V, 
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and references therein). Fly's Eye claims a value for -Eankic ~ 10^^'^ eV, a 
factor of three smaller than the AG AS A value (see Table V of [26]). This 
value is difficult to accommodate with the Z-burst mechanism, although one 
could perhaps gain agreement by increasing Zmax of the astrophysical neutrino 
sources. This would allow for more red-shifting of the nucleons, thereby low- 
ering the soft end of the Z-burst spectrum. A higher z^ax may be expected, 
for example, if the source is decaying massive particles or topological defects. 

The fit of the AGASA data with our total flux provides the normahzation 
of the photon and nucleon differential fluxes F, denoted as d(f)/dE in Figs. 2 
and 3 (in this case Fagasa = 6.6 x 10^"*^^ (m^ sr s)~^ -^pythia)- This allows 
us to determine the (assumed homogeneous and isotropic) flux of extreme- 
energy neutrinos close to the Z-resonance energy (more precisely, within the 
energy interval from -Bres/(1 + z:max) = ^res/3 to E^^es) to be 

F,,, ^ 6.6 X 10-=^«--4 . (3) 

eV m sr s 

This flux is shown in Figs. 2 and 3 , with the label 'Vee" • Within the level of 
accuracy of our simulation, we only determine the order of magnitude of this 
flux. Nevertheless, the flux requirement is formidable. This very-high energy 
flux is at present marginally ruled-out by present GLUE [30] and FORTE [31] 
bounds, as can be seen in Figs. 2 and 3 (see also [32]) The large required 
flux in the Z-burst hypothesis is also a challenge to theory. It appears that 
new physics is required to produce such a large neutrino flux at such extreme 
energy [33]. 

In this work, we have not included any effects of extragalactic magnetic 
flelds. For the analysis of this paper to hold as is, these fields should be 
sufficiently small, < 10~^ G. A way to weaken the bound derived here would 
be to assume the existence of large enough extragalactic magnetic fields, say 
10~^ G or larger, so that through electron pair-production and subsequent 
synchrotron radiation, photons would be eliminated as EECR primaries. 
This would reduce the event count above the AGASA end-point, and so allow 
a smaller neutrino mass. However in this case the flux of EE neutrinos at the 
Z- resonance energy, -Fl^^^, required to produce the observed EECR's would 
be larger. Moreover, the resonant energy, E^es oc l/rrii,, would be larger. 
Both the increased flux and the increased energy bring this possibility into 
conflict with the GLUE and FORTE bounds, on EE neutrino fluxes, shown 
in Figs. 2 and 3. Speciflcally, we flnd that without the photon contribution 
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to the Z-burst flux, we would require = 7.23 x 10 (eV sr s) ^ at 
Eres = 2.08 X 10^2 eV for = 0.2 eV, or F^^^ = 4.26 x 10-^6 (eV sr s)-^ 
at Ej-es — 4.16 X 10^2 eV for m,^ = 0.1 eV. Such possibilities are rejected by 
the GLUE and FORTE bounds. 

Another prediction of the Z-burst mechanism is the enhancement of the 
event rate at the GZK cutoff energy, Eqzk ~ 5 x 10^^ cV, due to the ac- 
cumulation of nucleons after energy-attenuation. Such an enhancement (or 
"GZK bump") is expected from any model which provides a nucleon spec- 
trum flatter than E'"^ above Eqzk (spectra falling faster than E~'^ do not 
provide sufficient flux above Eqzk to make an observable pile- up.). A hint 
for such a bump is seen in the AGASA data. Fig. 2 and 3. However, one 
must be careful in interpreting the shape as an enhancement. An alterna- 
tive explanation of the shape feature has been provided in [34] , namely that 
another energy-loss mechanism, p + j p + e'^ + e" , depletes the flux at 
energies just below Eqzk- 

Finally, a reliable prediction of the model is that most primaries above the 
ankle should be nucleons up to lO^^" eV or more, and photons and protons at 
higher energies. A much larger event sample is needed to test this prediction. 

5 Conclusions 

The non-observation of cosmic ray events at energies higher than ~ 2 x 10^° 
eV by AGASA provides a lower bound of about 0.3 eV on relic neutrino 
masses, if Z-bursts are responsible for the EECRs. This bound together with 
comparisons of the microwave background anisotropy to today's large-scale 
structure, leave a small interval for the relic neutrino mass, centered around 
0.3 eV, if Z-bursts explain the EECRs. The ~ 0.3 eV mass provides a very 
accessible signal for neutrinoless double beta-decay experiments (assuming 
neutrinos are Majorana particles), and is in fact compatible with a reported 
positive signal [35]. This ~ 0.3 eV mass interval may also be probed by the 
future KATRIN tritium end-point experiment [36] , and by expected improve- 
ments to the astrophysical bound from weak lensing of distant galaxies by 
intervening matter distributions [37], and from mining future Planck mea- 
surements of the CMB [38] . We note that our small allowed range for rrii, is 
incompatible with an upper bound of 0.15 eV that arises from invoking the 
see-saw mass spectrum as a source of successful early-Universe baryogene- 
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sis [39]. 

We have not included in this work any possible effects of cxtragalactic 
magnetic fields. For the analysis of this paper to hold as is, these fields 
should be sufficiently small, ^ 10~^ G. A survey of the literature shows 
that extragalactic fields much larger than 10~^ G are not expected, except 
in regions around galaxy clusters. In any event, large cxtragalactic magnetic 
fields which would reduce or eliminate the EE photon component of the Z- 
bursts (thus lowering the neutrino mass we obtain) would also lead to larger 
required EE neutrino fluxes at the Z-resonance energy. This possibility is 
rejected by the GLUE and FORTE bound. 

We have shown that Z-bursts may account not only for the super GZK 
events, but for the "ankle" and all EECR events above it. In doing so, 
the Z-burst mechanism meets the "naturalness" requirements of isotropy 
and spectral smoothness. In our simulation we found the "ankle" close to 
Eankie = 10^^-° cV as mcasurcd by AGASA. We found that most EECR 
primaries above the ankle should be nucleons up to about 10^*^ ° eV and nu- 
cleons and photons at higher energies. We also found that primary nucleons 
do accumulate at the GZK cutoff energy, which could account for the slight 
accumulation seen in the data (see, however, [34]). 

With rriiy ~ 0.3 eV needed to accommodate the AGASA end-point, the Z- 
burst hypothesis predicts a new absolute energy cutoff for EECRs, at E^cs — 
1.3 X 10^^ eV. However, due to the large multiplicity in Z-decay, the vast 
majority of incident primaries are expected to carry only a few percent of 
this cutoff energy. 

In summary, we believe that the Z-burst mechanism provides a plausible 
explanation to the puzzle of extreme-energy cosmic rays, not only for the 
super GZK events, but for the "ankle" and all events above it. But the 
hypothesis will stand or fall on the required nearness of nii, to 0.3 eV, and/or 
on the proximity of the required neutrino flux at E^es = 1.3 x 10^^ eV to 
existing flux limits. 

If Z-bursts are invalidated as the source of the observed EECRs, then we 
would encourage future, larger experiments (e.g.. Auger, EUSO) to continue 
the search for Z-burst signatures, for they provide the only known window 
to discovery of the rehc neutrino background. The Z-burst mechanism is 
entirely Standard Model particle physics coupled with Standard Hot Big- 
Bang Cosmology. Consequently, the mechanism itself cannot be ruled out. 
However, the Z-burst rate depends linearly on the resonant-energy neutrino 
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flux which Nature provides. Nature may not be generous. 
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